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ABSTRACT 

Elliptical, lenticular, and early-type spiral galaxies show a remarkably tight power- 
law correlation between the mass M, of their central supermassive black hole (SMBH) 



and the number Nqc of globular clusters: M, 



m. 



X N. 



1.08±0.04 

GC 



with m. 



1.7 X 



10 M0. Thus, to a good approximation the SMBH mass is the same as the total mass 
of the globular clusters. Based on a limited sample of 13 galaxies, this relation appears 
to be a better predictor of SMBH mass (rms scatter 0.2 dex) than the M,~a relation 
between SMBH mass and velocity dispersion a. The small scatter reflects the fact that 
galaxies with high globular cluster specific frequency Sn tend to harbor SMBHs that 
are more massive than expected from the M,-a relation. 

Subject headings: black hole physics - globular clusters: general - galaxies: star clusters 
- galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: evolution 



1. Introduction 



Supermassive black holes (SMBHs) have been detected in the centers of many nearby galaxies 



([Kormendv &: Richstondll995l : iMagorrian et al.lll998l : iGiiltekin et a 



correlated with several properties o f their host galaxies (INovak et al 



dispersion (the M.-g relat ion, e.g.. lFerrarese &: Merritt 



2001 



2002 



20091 ). The SMBH masses are 



20061') . in particular the velocity 



Gebhardt et al. 



20001 : 



Tremaine et al 



Giiltekin et al.ll2009l ) and the mass and luminosity of the spheroidal component — the entire 



galaxy in the case of ellipticals or the bulge in the case of lenticular and spiral g al axies ( Kormendv! 



1993 



Kormendv &: Richstone 



1995 



Marconi Hunt 



2003 



iticular and spiral g al axies (iKormendv^ 
Haring Rixl bood ) . ISpitler fc ForbesI 



((20091) find a tight correlation between SMBH and dark matter halo masses. As the dark halo 
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properties are inferred from the number of globular clusters in a galaxy this also indicates a con- 
nection between globular clusters and SMBHs. These correlations suggest a strong link between 
SMBH formation and galaxy formation, although the nature of this link is poorly understood. 

Numerous authors have investigated the possibility that the growth of galaxies and their 



SMBHs is regulated by their interactions (e .g., iHaehnelt fc Kauffmannll200d : iBurkert Sz Silkll2001 



Somerville et al.ll2008l : ICattaneo et al.ll2009l ). SMBHs grow by several mechanisms, including accre- 
tion of gas, swallowing stars whole, o r merging with other SMBHs acquir ed through a merger of their 
host galaxies. The Soltan argument (|Soltanlll982l : lYu Sz Tremainell2002l ) sugges ts that gas accretion 
is the dominant c ontributor to the SM BH mass budget, and both observations (jSanders et al.lll988l ) 
and simulations ([Hopkins et al.ll2006l ) suggest that much or most of this accretion occurs during 
mergers. SMBH growth through gas accretion can release substantial amounts of energy that can 
heat the interstellar gas, quench star formation, or even drive a wind that sweeps the galaxy free of 
gas, thereby halting star formation completely. Simulations of gas-rich galaxy mergers, including 
seed black holes, can reproduce the observed M,-a relation remarkably well given the simplicity 



of the empirical prescriptions used to model the accretion and feedback processes (ISpringel et al 



2005 



Cox et al 



2006 



Croton et al.l 120061 : iHopkins et al.ll2008l : [Johansson et al 



2009a 



The origin of the seeds of SMBHs is controversial. According to one hypothesis, the seeds were 
remnants of the first generation of metal-free, massive stars that formed at high redshifts. How- 
ever, the existence of bright quasars at redshifts of z > 6 demonstrates that SMBHs with masses 
exceeding 10^ Mq were already in place less than a billion years after the Big Bang. It is difficult for 
black-hole remnan ts from first-generation stars to grow fast enough to explain these observations 
(jMaver et al.ll2009l ). A second hypothesis is that the seeds are much larger ( "intermediate- mass" ) 
black holes of 10^-10^ M©, perhaps formed by the direct collapse of gas at the centers of proto- 
galaxies. 

Globular clusters (GCs) are among the oldest stellar systems in the universe and may have 
formed at the same time as the first stars. Their high stellar de nsities, somet i mes exceeding 
IO^Mm pc~^, lead to a variety of complex dynamical phenomena (jSpitzed 119871 : iHeggie &: Hut 



2003 



Binney Tremaind l2008l ) . Among these is mass segregation, through which heavy, com- 
pact, stellar remnants — neutron stars and black holes — spiral into the center by dynamical fric- 
tion. Once these arrive at the center i t is possible, though far from certa i n, that they merge to 



form an intermediate-mass black hole (Lee 



2004 



19871 : iQuinlan &: Shapiro 



1987 



Portegies Zwart et al 



Kawakatu &: Umemurall2005l ). There is significant observational evidei ice for intermediate- 



mass BHs with masses 4 x 10^-4 x IO^A^q in the centers of several GCs (IGerssen et al 



2002 



Gebhardt et al.l 1 20051 : iNoyola et al.l l2008l : Ivan der Marel AndersonI l2010l ) , but this evidence is 
still controversiao 



^ An additional uncertainty is whether some of these systems might be tidally stripped dwarf galaxies masquerading 
as GCs. 
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The number of globular clusters in a galaxy, Nqc^ is roughly proporti onal to the total luminos- 
i ty of the galaxy's spheroidal component. This relation was quantified by [Harris &: van den Bergh 
(|l98ll ). who introduced the specific globular cluster frequency Sn-, defined as the number of GCs 
per unit absolute visual magnitude My = — 15, 



Sj, ^ Ngc X l00-4(*^^'+i5) 
where My is the magnitude of the spheroidal component. 



(1) 



Brodie &: Stradeij (|2006l ) have summarized the progress that has been made in the quarter- 
century since the work by Harris & van den Bergh. It has become clear that star cluster populations 
are powerful tracers of galaxy evolution and that the observed correlations between globular clus- 
ter and galaxy properties provide valuable information about thei r joint formation. One of the 
most comprehensive studies of early- type galaxies is bv iPeng et al.l (j2008l ). who measured specific 
frequencies for the globular cluster systems of 100 elliptical and lenticular galaxies in the Virgo 
cluster. They find that early-type galaxies with intermediate luminosities (—22 < My < —18) 
typically have Sn ~ 1.5, while lu minous galax ies have Si\f ~ 2-5. The dominant galaxy M87 has 
an even larger specific frequency (jRacing |l968|) , estimated by Peng et al. to be S'at ~ 13. 



The formation of GCs is not well- understood (see, e.g., iBrodie &: Straded l2006l for a review). 
An important clue is that gas-rich merging gal axies contain large numbers of young massive star 
clusters that presumably formed in the merger ( Schweizer 1987 : Whitmore &: Schweizerl 19951). As 
this p opulation of clusters ages it is likely to evolve into a population of "normal" GCs (jFall &: Zhang 
200 ll ). Another scenario is the combined formation of SMBH seeds and globular clusters in super 
star- forming clumps of gas-rich galactic disks at z ~ 2 ([Shapiro et al.ll20ld : iMcLaughlin &: Pudritz 



19961 ). 



In summary, (i) both the SMBH mass M, and the total number of GCs Nqc are roughly 
proportional to the total luminosity of the spheroidal component in early-type galaxies; (ii) GCs 
may provide the black-hole seeds from which SMBHs grow; (iii) both the growth of SMBHs and the 
formation of GCs appear to be associated with major mergers or global gravitational instabilities 
in gas-rich protogalaxies. Given these observations, it is natural to ask how the properties of the 
GC population in early-type galaxies are correlated with the properties of their associated SMBHs. 

In this paper we show that there is a tight, power-law relation between the mass of SMBHs 
and the total number of globular clusters in elliptical, lenticular and early-type spiral galaxies. 
Remarkably, this relation appears to have even less scatter than the classic relation between SMBH 
mass and the velocity dispersion of the host galaxy. The relation can be approximately characterized 
by the statement that the SMBH mass equals the total mass in GCs. 



- 4 - 



The BH-GC relation 



For this analysis we have selected all elliptical, lenticular and early-type spiral galaxies with 
good estimates of the SMBH mass M, and the total number of globular clusters Nqc- 

(I2OO9I), 



Most of the SMBH masses were taken from the recent compilati on by 



Giiltekin et al 



wit h revised values for NGC 4486 (M87) and NGC 4649 (M60) from iGebhardt fc Thomas! (|2009, ) 
and lShen fc Gebhardtl (j20ld ). respectively. We also include two unpublished SMBH measurements, 
provided by Karl Gebhardt: NGC 4472 (Shen, J., Gebhardt, K. et al. 2010, in preparation) and 
NGC 4594 (Gebhardt, K., Jardel, J. et al. 2010, in preparat ion). Most of the GC numbers were 
taken from the AC S Virgo Cluster Survey (|Peng et al.l |2008| . first choice) or the compilation by 
Spitler et al.l (j2008l . second choice) . Table 1 summarizes the data; the Table contains entries for the 
Milky Way, the Sb spiral M31, and the dwarf elliptical M32, but these are not included in the fits 
below. In three cases (NGC 1399, NGC 3379, and NGC 5128) where there are two good-quality 
values in the literature that differ by two standard deviations or more, we have included both 
estimates in our fits, each at half weight. Note that errors in GC numbers are significantly larger 
than Poisson, due to uncertainties in the extrapolation of the GC luminosity functions, background 
subtraction, and corrections of the observed area of the galaxy to the whole system; the last of 
these is a particular concern because the radial distribution of GCs does not always follow the 
radial distribution of light. 

The points in Figure [T] show SMBH mass M, as a function of GC population Nqc- We fit this 
data to an assumed underlying relation of the form logM, = a + f3log{NGc) (all logs in this paper 
are base 10). We determine the best-fit values o f a and by minimizin g including errors in both 
observational parameters, using the methods in lTremaine et al.l (|2002l ). There exists a surprisingly 
tight correlation (dashed line). 



log 



M, 







i.U ± 0.04) + (1.08 ± 0.04) log 



Ngc , 
500 ' 



(2) 



the psr degree of freedom is 6.6. For comparison, the M,-a relation for the same sample, shown 
in the upper left panel of Figure [21 is 



log 



M, 

Mo 



.36 ±0.04) + (4.57 ±0.25) log 



a 



200 km s" 



(3) 



with psr degree of freedom of 8.5. Thus, in this admittedly small sample {N = 13), the corre- 
lation of SMBH mass with globular-cluster number is actually tighter than the classic correlation 
with velocity dispersion. 

We have also carried out unweighted fits in which we ignore the observational errors in disper- 
sion, luminosity, and GC number and minimize the rms residual e in the log of the SMBH mass 
(weighted by the observational errors in mass). For the mass-dispersion relation, e = 0.30 dex and 
for the mass-luminosity relation e = 0.38 dex, while for the mass vs. GC number relation e = 0.21 
dex, substantially smaller than the other two. We used the same procedure to fit to a relation of 
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the form logM, = a + /3ilogcJ + f32^ogL + P^log Nqc and find that e = 0.19, only marginahy 
smaher than the rms residual of 0.21 to the fit involving only Nqc despite the presence of the two 
additional free parameters /3i and /32- 

The smaller rms deviation e seen in the correlation between M, and Nqc implies that this 
is not a "secondary" correlation due to a tight correlation between Nqc and bulge properties, 
combined with a bulge-SMBH correlation. Additional evidence comes from the upper right panel 
of Figure [21 showing the correlation of Nqc with bulge visual luminosity Ly which is less good (x^ 
per degree of freedom of 34.8). We have also checked the correlation of Ngg with lumin osity and 



dispersion using the much larger sample of 62 galaxies in the sample of iPeng et al.l (|2008l ) that also 
have dispersions in HyperLeda and found fits of similar quality, with per degree of freedom of 
23 and 27 respectively. 

Another look at the data is provided in Figure [3l which compares residuals from the best-fit 
M,~(T relation (eq. [3]) on the horizontal axis to residuals from the best-fit Nqc^L relation 

w = ^ + ^^-^^ ^ ip^ 

on the vertical axis. In general galaxies with positive (negative) residuals in one quantity have 
positive (negative) residuals in the other. 

It is helpful to look at a few individual galaxies. Compare the galaxies M87 (NGC 4486) and 
Fornax A (NGC 1316), which have similar luminosities (My = —22.7 and —22.8 respectively). 
M87 has an SMBH mass that is larger than expected for its dispersion by 0.20 dex, while Fornax 
A contains an unusually small SMBH, by —0.42 dex. M87 has a specific frequency Sn = 12.2 
that is a factor of 3 larger than the average value for our sample of Sn = 4.0, while Fornax A 
is characterized by an unusually small value, Sm = 0.9. A second striking example is NGC 821, 
which has the largest residual from the M,-a relation in our sample — logM, is 0.82 smaller than 
predicted. This galaxy also has a small specific frequency, Sn = 1.3, not far from the low value 
seen for NGC 1316. 

There is growing evidence that the masses of black holes in th e most luminous (core ) gala xies 



have been underestima ted in some cases: including a d ark halo (iGebhardt &: Thomas! l2009l ) or 



allowing for triaxiality (jvan den Bosch &: de Zeeuwll2010l ) both tend to increase black- hole masses 



by a factor of two or so in the few luminous core galaxies modeled so far. To test for the effect 
of these revisions, we have increased the black-hole masses in core galaxies, by a factor of four if 
the models accounted for neither a dark halo nor triaxiality, and by a factor of two if the models 
accounted for one of these two effects. This change slightly steepens the M^-Nqc relation — the 
best-fit slope in equation ([2]) increases from 1.08±0.04 to 1.17±0.04 — and reduces the per degree 
of freedom from 6.6 to 4.8. It also reduces the psr degree of freedom for the M,-a relation, 
from 8.5 to 5.0. With this revised mass scale the unweighted best-fit M,-Ngc and M,-a give rms 
residuals of 0.20 and 0.26 respectively. 

Assuming a mean globular-cluster mass mac = 2 x 10^ Mq, equation ^ can be re-written in 



- 6 - 



terms of the total mass of the globular-cluster system Mqq = Nqc mcc, 

M, Mrr 

log = (0-14 ± 0.04) + (1.08 ± 0.04) log ■ 



o 



lO^M 



o 



(5) 



Remarkably, to a good approximation the mass of the SMBH is the same as the total mass of the 
globular clusters. 



3. Local Group galaxies 

The compact elliptical galaxy M32 has a central SMBH with mass M, = 3 x 10^ (Table 1) 
but no known globular clusters. This exception to the correlation between M, and Nqq observed in 
more luminous galaxies probably arises because t he globular clus ters in M32 have spiraled into the 
center of the galaxy through dynamical friction (jTremaind 119761 ) . The charac teristic inspiral time 
for a n object of mass m at initial radius in a galaxy with dispersion a is (jBinnev &: Tremaine 
20081 . eq. 8.12) 



1.65 rfa 
In A Gm 



1.2 X 10^° yr: 



In A V0.5kpc J \2x 10^ M 



-1 



75 km s 



(6) 



where In A is the standard Coulomb logarithm, m is normali zed to the mean g lobular-cluster mass 
of 2 X 10^ Mq and the dispersion is normalized to M32's ( Choi et al.ll2002l ). According to this 
estimate GCs within 0.5 kpc would have spiraled into the center within a Hubble time. M32's 
effective radius is only 0.14-0.18 kpc (jChoi et al.ll2002l ). so it is not surprising that most or all of 
the GCs in this galaxy have disappeared. 

The nearby Sb spiral galaxy M31 has 450 it 100 GCs (Table 1). According to equation ([2]) this 
leads to a SMBH mass of (1.2 it 0.4) x 10^ M©, in excellent agreement with the observed value of 
lAt^l X 10® M0. Thus the M,-Ngc correlation may work well for spiral galaxies as late as Hubble 
type Sb. 

The situation is less clear for the Milky Way, which contains a small SMBH of mass (4.3lg 4) ^ 
10^ M0 (Table 1). According to equation ([2]) this galaxy should contain only about 20 GCs compared 
to the observed population of 160 zt 20. The Milky Way is quite different from the elliptical, 
lenticular, and early-type spiral galaxies that we have discuss ed so far: it has a significantly later 
Hubble type (Sbc) and its bulge is probably a pseudobulge (jBinneyl l2009l : IShen et al.l l20ld ). In 
fact, most of the GCs in the Milky Way are associated with the Galactic halo, not the bulge, and 
as SMBH mass appears to correlate with bulge luminosity rather than total luminosity in spirals, 
it might be reasonable to expec t that in late-type spirals equation ([2|) applies only to the number of 
bulge GCs. iForbes et al.l (|200ll ) argue that metal-rich GCs at radii < 5 kpc are mostly associated 
with the bulge, and estimate that there are 35 it 4 such clusters in the Milky Way. If we insert 
this number in equation ([2]) we obtain an SMBH mass of (7.8 ± 1.2) X 10^ M0; this is somewhat 
larger than the observed mass but in reasonable agreement since some of the metal-rich GCs in 
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this region are likely to be disk GCs. We conclude that the Milky Way offers no strong evidence 
for or against the proposed correlation. 



4. Discussion and conclusions 

We have found that there is a strong correlation between the number of globular clusters and 
the mass of the central SMBH in early-type galaxies. This correlation appears to be at least as tight 
as the well known correlation between velocity dispersion and SMBH mass, although this conclusion 
is based on only 13 galaxies. To a reasonably good approximation, the BH-GC correlation simply 
says that the mass of the central SMBH in an early-type galaxy is equal to the mass of its GCs 
(eq. [5]). We suspect that the proportionality of the SMBH mass to the total globular-cluster mass 
offers insight into their formation processes, but the near-equality of the masses is a coincidence. 

Most galaxies have GC populations with a bimodal color distribution: there are r ed (metal-rich) 



and b lue (metal-poor) peaks, presumably reflecting two sub-populations of GCs (e.g.. lBrodie fc Strader 



20061 ). It is interesting to investigate whether the SMBH mass is correlated with one or the other 



of these su b -popu latio ns. Table 1 shows the red cluster fraction /red for 11 galaxies, taken from 



Peng et al.l (|2008l ) and iRhode fc Zepa (|2004l ). Note that /red is rather constant, with mean and 
standard deviation 0.3 it 0.1. The lower right panel of Figure [2] shows the M,-Ngc correlation sep- 
arately for the blue (triangles) and red (circles) clusters. As expected from the small rms variation 
in the red cluster fraction both correlations are of similar quality. 

The origin of the M,-Ngc relation is obscure. One possibility is that both the growth of 
SMBHs and the formation of GCs are associated with major mergers, so that galaxies that ex- 
perienced a recent major merger will have anomalously large SMBH masses and GC populations. 
Another possibility is the correlated formation of SMBH seeds and globular clusters in gas-rich 
young galaxies. 

An important next step is to expand the sample of galaxies having both reliable SMBH masses 
and reliable globular-cluster populations. 
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Fig. 1. — The number of globular clusters I^gc is shown as a function of SMBH mass Mbh for the 
13 giant elliptical, lenticular and early- type spiral galaxies in Table 1. Open circles connected by 
dotted lines denote the galaxies NGC 1399, NGC 3379 and and NGC 5128 for which two estimates 
of the SMBH mass are given. The dashed curve shows the fit given by equation ([2]). The location 
of M31 is also plotted as an open triangle, but this galaxy does not contribute to the fit. 
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Fig. 2. — The upper and lower left panels show the correlation between M, or Nqc, respectively, 
and the velocity dispersion. The dashed curve in the upper left panel corresponds to equation 
dSD, and the line in the lower left panel is log(iVGc/500) = 0.15 + 4.541og(cr/200 km s~^). In the 
upper right panel Nqc is plotted versus the visual luminosity of the host galaxy; the line shows the 
correlation given by equation In the lower right panel blue triangles and red points correspond 
to the number of blue and red globular clusters, respectively, versus SMBH mass; in this panel the 
galaxies NGC 1399, NGC 3379 and NGC 5128 are represented by the geometric mean of the two 
SMBH mass estimates. 
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Fig. 3. — Residuals from the best-fit relation between SMBH mass and velocity dispersion (eq. [3]) 
and the best-fit relation between number of globular clusters and bulge luminosity (eq. U]). Dashed 
lines and open circles denote the galaxies NGC 1399, NGC 3379 and NGC 1399 for which two 
estimates of the SMBH mass are given. 
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Table 1: Properties of galaxies with known SMBH masses and globular-cluster populations. The columns give morphological type, SMBH 
mass Mbh, number of globular clusters Nqc, velocity dispersion cr, absolute visual magnitude My for the entire galaxy and for the 
spheroidal component alone as determined from B/T, globular cluster specific frequency Sm relative to the spheroidal component (eq. [TJ, 
red cluster fraction /red and bulge/total luminosity ratio B/T. We list t he Local Group galaxi es M 31, M32, and the Milk y Way below 
the ho rizon tal line, but these a re no t included in our fits. References: fl) iGiiltekin et al.l |2009l) : (2) IPeng et al] (|20Q8[ ): f3) ISpitler et~ar 



(|2008l ): (4) iGomez et al.l (|200l[ ): (5) iNowak et al.l (120081): f6a) Shen J., G ebh ardt, K. et al. 201 (to b e su bmitted); (6b) Gebha r dt, K ., 
Jar deh J. et al. 2010 (t o be submitted): (7) iGebhardt fc ThomasI (l2009h: (8 ) Ishen fc Ge bhardt! (l2010l): (9) IXundu k. Whitmorel (l200ll): 



Forbes et al 



mO^: (18) R hod e fc ZeDjI 120041: (14) Ivan den Bosch de 



Giiltekin et al.l (|2009f ) for NGC 4594, and lKormendv et al 



Zeeuw 



(201 



('in) [ Bender e t al.l hoO^: (11) iGillessen et al.l hm^-. (12) 

(|2010|). B /T ratios are ta ken from lKormendv et"air |2009l) for NGC 4564 , ^ ^ , 

for M31. ' jCiiltekin et al.l |2009t ) give two estimates for the SMBH masses in NGC 1399 and NGC 5128, and give half-weight to each in then 
fits. We adopt a similar procedure: we include both estimates but increase the error bars on each by ^2 so that their combined contribution 
to the fit is the same as that of a single galaxy. ''The Milky Way contains a pseudobulge. 



